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We present preliminary results on the proton spin structure function in full QCD. The measurement has been 
done using 4 flavours of staggered fermions and an improved definition of the lattice topological charge density. 



1. INTRODUCTION 

The so-called proton spin crisis is still calling 
for a theoretical explanation in the context of 
QCD. The problem may be stated as follows. Let 

be the singlet axial current 



Nf 



(1) 



Its on-shell nucleon matrix clement may be writ- 
ten as 

(jp,s\jl\p',s')^ 

u{p,s) [Gi(fc2)7^75 - G2{k')k^j,] u{p\s'), (2) 

where k = p' — p. In the naive parton model 
AS = Gi(0) can be interpreted as the fraction of 
the nucleon spin carried by the quarks, and so is 
expected to be close to 1. On the contrary, ex- 
perimental determinations [Q-^ lead to a value 
AS ~ 0.2, which is also in contrast with the OZI 
expectation AS ~ 0.7 Q]. It is the small exper- 
imental value of AS which is referred to as the 
proton spin crisis. 

A numerical simulation on the lattice is the 
ideal nonperturbative tool to obtain an estimate 
of AS from first principles. A direct computation 
of the matrix element of eq. (||) involves the eval- 
uation of a disconnected diagram, characterized 
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by a very noisy signal (5|-0. We shall adopt a 
different method P-[lO|] , which makes use of the 
axial anomaly 



d^^jlix) - -2NfQix) 

Nf is the number of flavours and 

^2 
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(4) 



is the topological charge density. The on-shell 
nucleon matrix element of Q{x) can be written, 
using Eq. (||), as 



{p,s\Q{x)\p',s') 



rriN 



f 



A{k'^)u{p,s)i-/5u{p',s'), 



(5) 



where toat is the nucleon mass and A{k'^) = 
Gi(fc2) + G2(fc2)fc2/mAr, whence A{0) = AS. 
The last relation does not hold in the quenched 
approximation, since G2{k^) develops a pole at 
P — in that case ||^ . 

Therefore, AS can be calculated on the lat- 
tice by evaluating the nucleon matrix element of 
the topological charge density Q{x) in full QCD. 
However, a necessary condition to do this is to 
have the correct sampling of topological modes. 
Here we will give some preliminary results ob- 
tained using the standard HMC algorithm with 
four flavours of dynamical staggered flavours and 
will show that the feasibility of the determina- 
tion is affected by the bad sampling of topological 
modes obtained with the HMC algorithm [p| . 
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2. THE METHOD 

We have simulated four degenerate flavours of 
staggered fermions at /? = 5.350 on a 16"^ x 24 
lattice, using the standard HMC algorithm de- 
scribed in Ref. and the Wilson action for the 
pure gauge sector. A trajectory length r — 0.6 
was used, with time step = 0.004. The value 
a ■ ruf = 0.01 was chosen for the fcrmion mass. 

In order to evaluate the matrix element of Eq. 
(^) we need a lattice regularization Ql{x) of the 
topological charge density. The renormalization 
properties of Ql{x) are nontrivial [|l3|: it mixes 
with the continuum Q{x), as well as 

with p{x) — 'J2fJi "^iV'iTsV'i- However, it can be 
shown that the last two operators can be safely 
neglected in a first approximation [|3|, so that 
the following relation holds 

Ql =a%Q + 0(a^). (6) 

We use the improved smeared operators Q l {x) 
proposed in Ref. and in particular the op- 

erator defined at the second smearing level. The 
value of Zq can be obtained non-perturbatively 
using the so-called heating method |]l5|,^: the 
short range fluctuations responsible for the renor- 
malization are thermalized by applying a few up- 
dating steps on a discretized classical configura- 
tion of charge 1. The thermalization shows up as 
a plateau in {Ql) plotted against the number of 
heating steps performed, as is shown in Fig. 1. 
Zq is then measured on this plateau. We have 
obtained Zq = 0.56 ± 0.04 with the values of the 
parameters used in our simulation. 

The next step is to compute the correlation 
function 

= ^ OM \ I 0)P±Qi(p, r)i?(p, <)), (7) 

2NfZQ\p\ 

where 

QL(P,r) = ^e^P^Qi(a;,T), (8) 

X 

and P± is the spin projection operator. For B 
we have taken the wall source, after performing a 
Coulomb gauge fixing on the i = time slice by 
means of a steepest-descent algorithm. For B we 
have taken the usual local baryon operator | lOp^ 

B = eabcX'x'x'- (9) 
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Figure 1. Determination of Zq for the 2-smeared 
operator. The number of heating steps performed 
is measured in terms of molecular dynamics time. 



The following behaviour is expected for 

C{t) 00] 

C{t) = A+ AE e-™"^-^«(*-^) + 

A. ASa (-l)te~™'^^--^'^(*-^\ (10) 

where ^+(-) and TOAr(A) can be determined by a 
fit to the full baryon propagator {B{0)B{t)) ||l^, 
H^. AEa and rriA are the axial baryonic charge 
and mass of the A(1405) particle respectively [T7[ |. 

3. RESULTS AND CONCLUSIONS 

We have collected a sample of 230 configura- 
tions, each separated by 15 HMC trajectories, us- 
ing 700 hours of a 25 Gflops APE/QUADRICS 
machine. 

In Fig. 2 we show the topological charge distri- 
bution in our sample. It is clear that the correct 
thermalization for the topological charge has not 
yet been achieved: (Q) ^ and the distribution 
is not symmetric under Q —Q. The time his- 
tory of the topological charge is shown in Fig. 3, 
where the inefhciency of the HMC algorithm to 
change the topological sector is clearly visible. 

Despite this deficiency in our ensemble of con- 
figurations, we have attempted an estimate of AS 
and ASa. From a fit to the full baryon propaga- 
tor we have obtained 

A+= 0.152(12), = 0.693(11) 

= -0.192(35), mA = 0.797(25) (11) 
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Figure 2. Topological charge distribution after 30 
cooling steps. 
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Figure 3. Time history, in units of molecular 
dynamics time t, of the topological charge Q. 



We have thereafter measured the correlation func- 
tion of Eq. (0). The zero momentum transfer 
limit has been approximated by using |p| — 27r/16 
(~ 600 MeV on our lattice). An average has been 
performed over the three possible spatial direc- 
tions for p and over the two possible polariza- 
tions. 

We have observed a poor signal for C{t) and, 
by combining data at t — A and r — 20, we have 
obtained 

AE- 0.04 ±0.04, AEa ~ 0.05 ±0.05 (12) 

The error analysis has been done by a jack- 
knife method. However, the error quoted in Eq. 
( p!2| ) does not take account of the systematic er- 
ror coming from the bad sampling of topologi- 
cal modes. By looking at the effective number 
of topological sectors explored during the simu- 
lation (Fig. 3), we roughly estimate the effective 
error to be up to 4-5 times the one quoted in Eq. 

Therefore our data can only be regarded as pre- 
liminary and no conclusive results can be pro- 
duced unless a correct sampling of topological 
modes is obtained. This goal seems to be hardly 
feasible for the standard HMC algorithm: the 
possibility of achieving it by using different al- 
gorithms is currently under investigation. 
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